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Abstract. Recent data from ATIC, CREAM and PAMELA revealed that the energy spectra of cosmic 
ray (CR) nuclei above 100 GeV/nucleon experience a remarkable hardening with increasing energy. This 
eifect cannot be recovered by the conventional descriptions of CR acceleration and diffusive propagation 
processes. Using analytical calculations, I show that the hardening effect can be consequence of a spatial 
change of the CR diffusion properties in different regions of the Galaxy. I discuss the implications of this 
scenario for the main CR observables and its connections with the open issues of the CR physics. 



1. Introduction 

Understanding the origin of the cosmic ray (CR) energy spectrum is central to astrophysics and has long 
been the focus of intensive study. The spectrum of CR nuclei at ~ 10^-10^ GeV/nucleon is thought to be 
the result of diffusive shock acceleration (DSA) mechanisms in supernova remnants (SNR), followed by 
diffusive propagation in the interstellar medium (ISM) [1]. The conventional descriptions predict source 
spectra such as E~''' for primary nuclei {e.g., H, He) which are steepened as E^^^^ by diffusion. The 
data constrain v + 5 ^2.1 (depending on the element) and 5 ~ 0.2-0.7, whereby v ~ 2.0-2.5. In DSA 
calculations, the source spectral slope should be v <2.2, which implies 5 >0.5. On the other hand, large 
values for 5 are at odds with the anisotropy observations at TeV energies, which favor 5 <0.3. On top 
of that, recent experiments ATIC-2, CREAM, and RA^MELA, reported a remarkable spectral hardening 
at energies above ~ 100 GeV/nucleon which cannot be explained by conventional mechanisms ||2j|3l|4l. 
Proposed explanations of this effect deal with acceleration mechanisms [5,6], nearby SNRs f7^ or multi- 
source populations lUllU. I propose here that the hardening originates from a spatial change of the CR 
diffusion properties in the different regions of the Galactic halo. The key hypotesis is that the diffusion 
coefficient K is not separable into energy and space terms as usually assumed. As I will show, this 
hypotesis leads to a pronounced change in slope for the energy spectra of CR nuclei at Earth, giving a 
good description of the present data. Remarkably, the model proposed here has also positive impacts on 
several open problems in the CR acceleration/propagation physics. 

2. The two-halo model 



I use a one-dimensional inhomogeneus diffusion model for CR transport and interactions [TOl. The 
Galaxy is modeled to be a disc (half-thickness h =100pc) containing the gas (number density 
n = 1 cm^"^) and the CR sources. The disc is surrounded by a diffusive halo (half-thickness L = 5 kpc) 
with zero matter density. For each stable CR nucleus, the transport equation is: 




(1) 



where N{z) is its number density as function of the z-coordinate, K{z) is the position-dependent 
diffusion coefficient and F'"*^' = (3cna^^'^^ is the destruction rate in the ISM at velocity /3c and cross 
section cj'"'^'. The source term Q is split into a primary term Qpn, from SNRs, and a secondary production 
term Qscc = X^j T^^^^^A'j, from spallation of heavier (j) nuclei with rate r^.P^"_ xhe quantities N, K, 
Q and F™''^ depend on energy too. Equation[T] can be solved in steady-state conditions {dN/dt = 0) 
with the boundary conditions A^(ibL) = ifTOl . For each CR nucleus, the differential energy spectrum 
is given by: 
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where Kq = K{z=0), A = Kq Jjj^' j^^, and A = X{z=L). In CR propagation studies, the diffusion 
coefficient is usually assumed to be separable as K{z,E) = f{z)KQ{E) (in homogeneus models 
f{z) = 1, whereby A = L). In these models, the functions A and A are independent on energy, so that the 
predictions at Earth (z=0) are spectrally uninfluenced by the choice of f{z). Conversely, I employ a non- 
separable K{z, E), which reflects the existence of different CR transport properties in the propagation 
volume, depending on the nature and scale distribution of the magnetic-field irregularities. In fact, while 
SNR explosions may generate large irregularities in the region near the Galactic plane, the situation in 
the outer halo is different because that medium is undisturbed by SNRs. From these considerations, the 
authors of Ref. [11] found that the turbulence spectrum in the halo should be flatter than that in the 
Galactic plane. This implies a strong latitudinal dependence for the parameter 5, which suggests spatial 
variations of the CR energy spectra. Noticeably, new data reported by Fermi/LAI on the diffuse 7-ray 
emission at ~ 10-100 GeV of energy seem to support these suggestions: the 7-ray spectra observed near 
the Galactic plane (latitude |6| <8°) are found to be harder than those at higher latitudes [ 12|. Following 
the above ai^guments, I adopt a simple two-halo model consisting in two diffusive zones. The inner halo 
is taken to surround the disk for a typical size of a few hundred pc (^ ~ 0. 1). Its medium properties 
are influenced by SNRs, which produce a steep turbulence spectrum in terms of energy density per wave 
number, w{K)dK ~ K^^^^dn, presumably close to the Kolmogorov regime 6 1/3. The outer halo 
represents a wider region, < \z\ < L, which is undisturbed by SNRs. Its turbulence spectrum is 
driven by CRs and should be flatter. For instance, CRs with rigidity spectrum R^'^ (where R = p/Z) 
excite turbulent modes of wave-numbers k ex. 1/R, giving a spectrum ^^k^^^^ , where 5' — 1 is of 
the order of the unity. This situation can be realized by a rigidity dependent diffusion coefficient of the 
type 

^) ^ ihHR/Ro)^ for\z\ < iL {inner halo) 
' {R/ Rq)^^^ for \z\ > {outer halo) , 

where ko = 0.04kpc/Myr specifies its normalization at the reference rigidity Rq =5GV. For R > Rq, 
the diffusion coefficient of Eq.[3] produces a higher CR confinement in the inner halo (with 6 ~ 1/3), 
whereas the outer halo (with A = 6' — 6 ~0.5 - 1) represents a reservoir from which CRs leak out 
rapidly and can re-enter the inner halo. The non-separability of K has a remarkable consequence on the 
model predictions at z = 0. It can be understood if one neglects the term F™'^' and takes a source term 
Qpri oc R^'^. From Eq.|2]and using Eq.[3l one finds 

Jo = j(z = 0) ~ ^ {e {R/Ro)-"-' + {l-0 {R/Ro)-"-'-^} , (4) 

which describes the CR spectrum as a result of two components. Its differential log-slope as a function 
of rigidity reads 

91ogJo A 

l{R) = — tt; ^ u + -\ T T- , (5) 

d^ogR i + _L.{R/R^f' 




Figure 1. Left: CR spectra of H and He. Right: mean logarithmic mass and B/C ratio. Model 
calculations are compared with the data E[2llll3llllll5l|l6l|l7l[l8l[l9l|20l|2ll. 



which shows a clear transition between two regimes. In practice the low-energy regime (7 « + (5 + A) 
is never reached due to spallation, that becomes relevant and even dominant over escape (r™*^' > ^). In 
this case the log-slope is better approximated by7~z^ + ^((5 + A) |[22l . The hard high-energy regime 
(7 « + (5) is determined by the diffusion properties of the inner halo only. In this limit one has A « ^L. 
The effect vanishes at all rigidities when passing to the homogeneus limit of ^ 1 (one-halo) or A ^ 
(identical halos), where one recovers the usual relation = v + S. From Eq.|2j it can be seen that the 
intensity of the harder component diminishes gradually with increasing \z\, i.e., the CR spectra at high 
energies are steeper in the outer halo, as suggested by the new Fermi/LAI data. 

Figure[T] shows the H and He spectra at Earth, from Eq.|2j at kinetic energies above lOGeV. The 
implementation of the model follows Refs. lITOl l23l . The source spectra are taken as Q oc R^"^, with 
u =2.3 for H and ly =2.16 for He El. The two halos are defined by L =5kpc and =0.5kpc. 
A Kolmogorov-type diffusion is adopted for the inner halo {6 = 1/3), while A is taken as 0.55, 
consistently with Ref. |11|. Results from the two-halo model (solid lines) are compared with those 
from the homogeneus model (dashed lines), which uses 6 =0.6, and with recent CR data. Below 
~ 10 GeV/nucleon, the solar modulation is apparent and it is described using a. force-field modulation 
potential =400MV 1.25,1 . Note that my model may be inadequate in this energy region due to 
approximations. Remarkably, my model reproduce well the observed changes in slope at > 100 GeV, in 
good agreement with the trend indicated by the data. It should be noted, however, that the sharp spectral 
structures suggested by the PAMELA data at ~ 300 GeV cannot be recovered. Figure[T]: shows the mean 
logarithmic mass, {ln{A)), which is described well by the two-halo model. Figure[T]l shows the B/C 
ratio from the two-halo model, which is predicted to harden as a K{E)/Kq{E), while its low-energy 
behavior is similar to that of the homogeneus model (with 5 = 0.6 in the whole halo). Interestingly, 
a slight flattening for the B/C ratio is also suggested by recent data, e.g., from TRACER. A multi- 
channel study of the AMS data will allow to resolutely test these features. The B/C ratio hardening is 
also connected with the global large-scale anisotropy amplitude, which should increases as oc Kq / A 
ll26l[T01 . Thus, the two-halo model predicts a gradual flattening of rj, as indicated by the present data. It 



is also interesting noticing that the anisotropy may be furtherly reduced at all energies if one accounts for 
a proper radial dependence for K ll27l . In summary, the model I proposed is potentially able to reconcile 
a weak energy dependent of the anisotropy amplitude (as suggested by observations) with relatively hard 
source spectra v ^ 2.2 (as preferred by the DSA theory), giving good fits to the B/C ratio. I recall that 
plain diffusion models employ u ss2.2 and 6 ssO.5 in the whole halo (which gives a too strong energy 
dependence for rj), whereas diffusive-reacceleration models use u ^ 2.4 and 6 ^ 1/3 (which are too 
challenging for the DSA and require strong reacceleration to match the B/C data). 

3. Conclusions 

I have shown that the spectral hardening observed in CRs may be consequence of a Galactic diffusion 
coefficient that is not separable in energy and space terms. From the model presented here, the hardening 
arises naturally as a local effect and vanishes gradually in the outer halo, where the CR spectra are also 
predicted to be steeper. This effect must be experienced by all CR nuclei as well as by secondary-to- 
primary ratios. Recent data from Fermi/LAI and TRACER seem to support this scenario, but the change 
in slope predicted by this model is more gradual than that suggested by PAMELA data. All these features 
can be tested with the data forthcoming by AMS, so that the present model will be discriminated against 
other interpretations. This scenario has an interesting impact on the CR physics. As shown, a Kol- 
mogorov diffusion for the inner halo (5 1/3) is consistent with relatively hard source spectra (u ~ 2.2), 
giving good descriptions of both the B/C ratio and the high-energy trend of the anisotropy amphtude. 
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